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We study the two orbital double-exchange model in two dimensions in the presence of antifer- 
romagnetic (AF) superexchange, strong Jahn- Teller coupling, and substitutional disorder. At hole 
doping X — 0.5 we explore the 'bicritical' regime where the energy of a ferromagnetic metal and 
a charge and orbital ordered (CO-00) CE state are closely balanced, and compare the impact of 
weak homogeneous disorder to that of a low density of strong scatterers. Even moderate homoge- 
neous disorder suppresses the CE-CO-00 phase and leads to a glass with nanoscale correlations. 
Dilute strong scatterers of comparable strength, however, convert the CE-CO-00 phase to a phase 
separated state with ferromagnetic and AF-CO-00 clusters. We provide the first spatial description 
of these phenomena and compare our results in detail to experiments on the half-doped manganites. 



The manganese oxides of the form Ai-a^A'a^MnOs in- 
volve a remarkable interplay of charge, spin, lattice, and 
orbital degrees of freedom [1 . This cross coupling is most 
striking in the half doped {x = 0.5) manganites many of 
which have a charge and orbital ordered insulating (CO- 
OO-I) ground state with 'CE' magnetic order - a zigzag 
pattern of ferromagnetic chains with antiferromagnetic 
(AF) coupling between them. The CE-CO-OO-I phase 
shows up in manganites with low mean cation radius 
{va) while systems with large va are ferromagnetic met- 
als (FM-M). The variation of va leads to a 'bicritical' 
phase diagram [2 with a first order boundary between 
the FM-M and the CE-CO-OO-I phases. 

Disorder has a remarkable effect on the bicriticality. 
Even moderate 'alloy' disorder, due to random location 
of A and A' ions at the rare earth site, converts the CO- 
00-CE phase to a short range correlated glass, but has 
only limited impact on the ferromagnet |2l |3l H]. The 
asymmetric suppression of spatial order by cation disor- 
der and the emergence of a charge-orbital-spin glass at 
low Va are one set of intriguing issues in these materi- 
als. Unusually, while alloy type randomness on the A 
site leads to a homogeneous glassy phase, the substitu- 
tion of a few percent of Mn (the 'B site') by Cr [H |6] 
leads to phase separation of the system [71 [H |9l [10] into 
FM-M and AF-CO-OO-I domains. The diflFerence be- 
tween A and B site disorder holds the key to the much 
discussed phase coexistence and spatial inhomogeneity in 
the manganites. 

In this paper we provide the first results on the relative 
effects of A and B type substitutional disorder on phase 
competetion in a manganite model. We study weak 'al- 
loy' disorder and dilute strongly repulsive scatterers. Our 
main results are: (i) Alloy disorder indeed leads to asym- 
metric suppression of long range order; moderate disor- 
der converts long range CE-CO-00 to an insulating glass 
with nanoscale inhomogeneities, while FM order is only 
weakened, (ii) A low density, ^ 4%, of strong scatterers 
in the CE phase leads to cluster coexistence of AF-CO- 
00 and FM regions and the ground state is a poor metal. 



(iii) The impact of strong scatterers depends crucially on 
whether they are attractive or repulsive, it correlates with 
the asymmetry of the 'clean' system about x = 0.5, and 
uncovers a new route for phase control. 

We consider a two band model for Cg electrons, Hunds 
coupled to t2g derived core spins, in a two dimensional 
square lattice. The electrons are also coupled to Jahn- 
Teller phonons, while the core spins have an AF superex- 
change coupling between them. These ingredients are all 
necessary to obtain a CE-CO-00 phase. We include the 
effect of disorder through an on site potential. 
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+ J^^^S,.S, -A^Q,.r, + ^^Q,2. (1) 
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Here, c and are annihilation and creation operators for 
Cg electrons and a, P are the two Mn-e^ orbitals dx2_y2 
and d^z^_r2, labelled (a) and (b) in what follows, are 
hopping amplitudes between nearest-neighbor sites with 
the symmetry dictated form: t^^ = t^^ = t, t^^ = t^^ = 

t-b = tba = -^/^/3, tl, = tl = t/V3, where x and y 
are spatial directions We consider effectively a lattice of 
Mn ions and treat the alloy disorder due to cationic sub- 
stitution as a random potential at the Mn site picked 
from the distribution PA{^i) = ^(^(^i — A) + ^(e^ + A)). 
The Cr doping case is modelled via Psi^i) = + 
(1 — r])5{ei), where r] is the percent substitution and V 
the effective potential at the impurity site. The Cg elec- 
tron spin is erf = 5];^^/ 4aa^aa'^io^(^'^ whcrc the F's are 
Pauli matrices. It is coupled to the t2g spin via the 
Hund's coupling Jh, and we assume Jn/t I. A is the 
coupling between the JT distortion = {Qix^Qiz) and 
the orbital pseudospin rf = Z^^^ 4a^r^^Q/3^, and K is 
the lattice stiffness. We set t = 1, K = 1, and treat 
the Qi and as classical variables [11 . The chemical 
potential fi is adjusted so that the electron density re- 
mains n = 1/2 which is also x = 1 — n = 1/2. For A 
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FIG. 1: Colour online: (a) The ground state at x = 0.5 for 
varying Jaf and A, in the absence of disorder, (b) The doping 
(n = 1 — x) dependence of the ground state for varying chem- 
ical potential /i and typical electronic couplings, A = 1.6 and 
Jaf = 0.1, near the FM-OD-CD k CE-CO-00 phase bound- 
ary. The phases in the vicinity of x = 0.5 are expected to show 
up in a cluster pattern on introducing disorder at x = 0.5. 



type disorder the mean value is = and the variance 
is A\ = ((e^ — e^)^) = A^, while for B type disorder 
e, = r]V and A^, = ((e, - e,)') = V^rji^ - r^). 

The clean CE ground state at x = 0.5 has been stud- 
ied earlier [12 , 13 , 14, 15 , 16 using mean field and Monte 
Carlo (MC) techniques and is well understood. The im- 
pact of disorder on the phase competetion appropriate to 
X = 0.5 has been studied on small clusters [iJl [181 Ull 
usually using simplified models either without orbital 
variables [18 or ignoring the electron-phonon coupling 
[19 . The difficulty of simulating the full model, eqn[l], 
on a large system has prevented any conclusive study. We 
use our recently developed travelling cluster approxima- 
tion (TCA) based MC [20 to solve the problem. Com- 
pared to exact diagonalisation (ED) based MC which can 
handle typical sizes ~ 8 x 8 we study the full model on 
lattices upto 40 x 40. In all our studies we use a moving 
cluster of size ~ 8 x 8 [20 to anneal the spin and phonon 
variables. 

Before discussing the effect of disorder we determine 
the clean ground state at x = 0.5 for varying Jaf and A, 
Fig.l.(a). 

At low A and low Jaf double exchange is the domi- 
nant interaction and kinetic energy optimisation leads to 
a homogeneous ferromagnetic state without any orbital 
or charge order (FM-OD-CD). This phase has a finite 
density of states at the Fermi level ep and is metallic. 
As Jaf is increased, keeping the JT coupling small, a 
magnetic state emerges with peaks in the structure fac- 
tor Smagisi) at q = {0, tt} or {7r,0} (we call this the A- 
2D phase), then an orbital ordered but uniform density 
CE phase, with simultaneous peaks at q = {0, tt}, {tt, 0}, 
and {7r/2,7r/2}. At even larger Jaf the dominant cor- 
relations are 'G type' with a peak at q = {7r,7r}. By 
contrast, increasing A at weak Jaf keeps the system fer- 
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FIG. 2: Colour online: (a) Experimental 'bicritical' phase 
diagram in the x = 0.5 manganites obtained for ordered and 
disordered (alloy) structures, (b) Our results: superposed 
phase diagrams at x — 0.5 for A = and A = 0.3. The long 
range CE-CO-00 for A > 1.55 at A = is completely wiped 
out at A = 0.3 while the FM-M phase at low A becomes an 
unsaturated FM with short range A- 2D type correlations. 



romagnetic but leads to charge and orbital order (FM- 
CO-00) for A ^ 1.6. Our interest is in a charge or- 
dered CE phase. Such a state shows up when both A 
and Jaf are moderately large. The TCA based phase 
diagram is broadly consistent with previous variational 
results [H [H [151 [E] and with ED-MC on small sys- 
tems [13]. 

Since the effect of disorder might be to create cluster 
coexistence [lU |22l of phases of different densities that 
arise in the clean limit, Fig.l.(b) shows the phases and 
phase separation windows that occur at a typical cou- 
pling, Jaf = 0.1 and A = 1.6. For these couplings the 
clean system is a CE-CO-00 phase at x = 0.5, a FM-M 
for X < 0.4, and an A-2D type AF for x > 0.55. 

To minimise the number of parameters, in what fol- 
lows we set Jaf = O-l- This is in the right ballpark 
considering the AF transition temperature at x = 1, and 
allows close proximity of the CE and FM-M phases. We 
mimic the bandwidth variation arising from changing rA 
by varying X/t across the boundary between CE-CO-00 
and FM-OD-CD, and now explore the effects of thermal 
fluctuation and disorder. 

The key experiment [2 on the effect of A site disorder 
on bicriticality compared an 'ordered' structure, where 
the rare earth and alkaline earth ions sit on alternate lay- 
ers^ with the 'disordered' case where they are randomly 
distributed. The result is reproduced in the left panel in 
Fig. 2. While the ordered case has large transition tem- 
peratures for the CO-00, CE, FM phases, e^c, a random 
distribution of A and A' ions destroy the CO-OO-CE 
phase and partially suppresses the ferromagnetic Tc. 

The right panel in Fig. 2 shows our result, where we 
superpose the clean phase diagram and the case with A 






FIG. 3: Colour online: Structure factors and resistivity at 
Jaf = 0.1 and A/t — 1.6. (a) The T dependence of the 
major peaks in the structure factor for spin order (SO), orbital 
order (00) and charge order (CO) in the clean limit. Note 
the clear separation of scales between Tco, Tsr and Tce- 

(b) The resistivity p{T) in the clean CE-CO-00 case and in 
the presence of A type and B type disorder, with Aa ~ As 
0.5. The As corresponds to V — 2 and dilution 77 = 0.08. 

(c) Variation of the major peaks in the magnetic structure 
factor with Aa at low temperature (T = 0.005). (d) Same as 
(c), now with B type disorder, V = 2 and varying 77. Note 
the emergence of the FM q = {0,0} peak around As = 0.4 
(ry = 0.04). 



type disorder = 0.3. In the clean limit at T = 
as \/t is increased there is a transition from a FM-M to 
the A-2D phase at \/t ~ 1.52, and then a transition to 
a CE-CO-00 phase at A/t > 1.55. On the FM-M side, 
\/t < 1.52, there is only a single thermal transition [23^ 
at Tc as one cools the system. At large A/t, however, 
cooling first leads to a CO-00 phase, at Tqoi without 
magnetic order, followed by strong features in Smag at 
q = {0,7r} and {0,7r}, showing up at Tsr^ indicative of 
stripelike correlations. Finally, at a lower T the system 
makes a transition to CE order. If we set t = 0.3eV, and 
use a factor of 3/2 to convert transition scales between 
2D and 3D, our Tc at bicriticality would be ~ 200K. 

In the presence of A type disorder with A^ = 0.3 we 
do not find any spatial order on the CE side in either 
the charge, or orbital, or magnetic sector, down to T ~ 
0.005. The absence of order in the CE-CO-00 side can 
be traced back to the 'random field' coupling directly to 
the charge order parameter n^. This breaks down charge 
correlations to the atomic scale. The ferromagnet being 
a q = state is more robust to A type disorder [18] . 

There are short range stripelike magnetic correlations 
that persist as peaks at q = {0, tt} and {tt, 0} in 5'ma^(q)- 
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FIG. 4: Colour online: MC snapshots and magnetic structure 
factor at low temperature, T — 0.01, size 40 x 40. Left row: 
A = 1.6, non disordered, middle row, A = 1.6, A type disor- 
der with Ae// = 0.5, right row, A = 1.6, B type disorder with 
y = 2, ?7 = 8%, Ae// = 0.56. Top panel shows the nearest 
neighbour magnetic correlation 8^.8^+5, where 5 — x ov y. 
Middle panel shows the charge density {rii) for the configura- 
tion above. Bottom panels shows the MC averaged Smag{^)- 
In each panel q = {0,0} at the bottom left corner, q = {tt, 0} 
at the bottom right corner, etc. 

The onset of this feature is shown by the (red) dotted line 
in Fig.2.(b). This appears even on the ferromagnetic side 
below The Tc itself is somewhat suppressed by dis- 
order and the ground state is an unsaturated ferromagnet 
(u-FM). Our analysis of the structure factor in the disor- 
dered system, however, does not suggest any coexistence 
of two distinct locally ordered phases at any A. A type 
disorder in the bicritical regime does not induce phase 
coexistence. We have confirmed this directly from the 
spatial snapshots as well, as we discuss later. 

We have explored A type disorder with strength A^ = 
0.1, 0.2, 0.3 and 0.4, over the range of X/t shown in 
Fig.2.(b). We now specialise to X/t = 1.6, which is a CE- 
CO-00 phase near the clean phase boundary in Fig.2.(b) 
and explore the impact of A type and B type disorder in 
detail. Fig. 3. (a) shows the T dependence of the major 
peaks in the spin, charge and orbital structure factor in 
the clean limit at X/t = 1.6, for reference, illustrating the 
distinct Tco^ Tsr and Tce scales. 

The naive expectation is that disorder would lead to 
cluster coexistence [211 E2] of AF-CO phases, that arise 
for X > 0.5, with the FM-M phase at x < 0.4, Fig.l.(b). 
Fig.3.(c) shows how the peaks in Smagi^) evolve with 
A A at low temperature (T = 0.01). The peak at 
q = {7r/2,7r/2} vanishes quickly, leading to a phase with 
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stripelike correlations, and the q = {0,7r}, {tt, 0} peaks 
also vanish for > 0.6 leaving a glass. The response 
to B type disorder is more interesting. We have explored 
y = 1, 2 and 4 and = 2, 4 and 8%. Since Cr is believed 
to be in a t^g^^ state we focus here on V = 2 which is 
sufficiently repulsive to force (rii) = (e^ state) at the 
impurity sites. The response, as we vary the fraction of 
scatterers (ry), is similar to A type at weak A^- However, 
before the peak at q = {0, tt}, {tt, 0} vanishes we see the 
emergence of a peak at the ferromagnetic wavevector, 
q = {0, 0}. There is a window at intermediate r] where B 
type disorder leads to coexistence of FM and CO-OO-AF 
regions. In terms of transport, Fig.3.(b), intermediate 
A type disorder strengthens the insulating character in 
p(T), while B type disorder of comparable variance leads 
to an insulator-metal transition on cooling, and a (poor) 
metallic state at low temperature. 

The top row in Fig. 4 compares low temperature MC 
snapshots of the magnetic correlations in the clean sys- 
tem at A = 1.6 (left), to that with A^ = 0.5 (center) and 
Ab = 0.56 (right). The respective panels in the middle 
row show the electron density (n^) corresponding to the 
panels above. The panels at the bottom are the thermally 
averaged Smag{^) in the three cases. In the clean limit 
the magnetic correlations are CE, with a checkerboard 
density distribution, and simultaneous magnetic peaks 
at q = {0,7r}, {tt, 0} and {7r/2,7r/2}. For A type disor- 
der there are stripelike magnetic correlations with small 
(atomic scale) FM clusters but no signature of phase co- 
existence. The density field is also inhomogeneous in 
the nanoscale, with only short range charge correlations, 
and Smagi^) has weak peaks at q = {0,7r} and {tt, 0} 
but no noticeable feature atq = {0,0}. B type disorder, 
however, leads to FM regions coexisting with stripelike 
AF correlations. The density field shows a correspond- 
ing variation, being roughly homogeneous within the FM 
droplets (with local density n ~ 0.6), and a CO pattern 
away from the FM regions. Smag{^) now has peaks at 
q = {0, tt}, {tt, 0} and {0, 0}, as seen earlier in Fig.3.(d). 

We explain the difference between the impact of A type 
and B type disorder as follows. (1) The introduction 
of A type disorder does not lead to coexistence of large 
FM-M and AF-CO-00 clusters, despite the presence of 
a PS window in the clean problem, Fig.l.(b), because 
(a) atomic scale potential fluctuations disallow CO co- 
herence beyond a few lattice spacings, while (b) homo- 
geneous FM-M clusters are destabilised by the disorder 
and become charge modulated. The result is a nanoscale 
correlated insulating glassy phase. (2) Dilute strongly 
repulsive scatterers act very differently: (a) they force 
an state at the impurity sites and generate an 'excess 
density' 0.5 x r] which has to be distributed among the 
remaining Mn sites, (b) the parent x = 0.5 CO phase can- 
not accommodate this excess charge homogeneously and 
the system prefers to phase separate into x ~ 0.5 AF-CO 
and X ^ 0.4 FM clusters, (c) unlike the A type case, the 



FM clusters can survive and percolate since at low r] there 
can be large connected patches without a B type site. We 
have verified this explicitly for several impurity configu- 
rations. Making the B site potential strongly attractive 
leads to a glassy AF-CO state since carrier trapping re- 
duces the effective electron count and forces the system 
towards a combination of x > 0.5 phases in Figl.(b). 

In conclusion, we have reproduced all the key effects of 
A and B type disorder on phase competetion in the half 
doped manganites. Our results suggest that B site im- 
purities can be chosen to engineer phase control and the 
percolative conduction paths can be controlled through 
choice of dopant locations. 
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